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Abstract. We study the Aharonov-Bohm-Coulomb problem in a graphene ring. 
We investigate, in particular, the effects of a Coulomb type potential of the form 
£/r on the energy spectrum of Dirac electrons in the graphene ring in two different 
ways: one for the scalar coupling and the other for the vector coupling. It is 
found that, since the potential in the scalar coupling breaks the time-reversal 
symmetry between the two valleys as well as the effective time-reversal symmetry 
in a single valley, the energy spectrum of one valley is separated from that of 
the other valley, demonstrating a valley polarization. In the vector coupling, 
however, the potential does not break either of the two symmetries and its effect 
appears only as an additive constant to the spectrum of Aharonov-Bohm potential. 
The corresponding persistent currents, the observable quantities of the symmetry- 
breaking energy spectra, are shown to be asymmetric about zero magnetic flux in 
the scalar coupling, while symmetric in the vector coupling. 
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1. Introduction 

Graphene is a two-dimensional composite lattice with honeycomb structure and 
consists of two sublattices of carbon atoms. The corresponding Brillouin zone (BZ) 
is also a hexagon with high symmetry points at vertices as well as center and side 
(see Fig. [I]) [TJ [2]. Near the vertices of the BZ the low-energy electronic spectra are 
linearly dependent on the magnitude of momentum, forming conical valleys, and the 
dynamics of electrons in graphene can be formulated by Dirac equation of massless 
fermion [2111]. Because of this Dirac fermion-like behavior of electrons, graphene may 
offer an opportunity to test the various predictions of the planar field theories [SI [5] by 
experiment with solid state material. In this respect there have been much effort to 
connect the two-dimensional (2D) field theory with the graphene physics [T;, including 
the Aharonov-Bohm (AB) effect [8, in a ring geometry of graphene [9| flO l ITT | [12], 

The studies of the AB effect in a graphene ring were mostly concerned about 
the breaking of time-reversal symmetry (TRS) that yields a splitting of the valley 
degeneracy due to the time-reversal symmetry in graphene. In this regard, the authors 
of Ref. [5] have demonstrated that the splitting of the degeneracy in a single valley can 
be controlled by the threaded magnetic flux and a confinement potential of the Dirac 
electrons on the AB ring of graphene. An interesting result from the work was that 
the confinement potential, introduced as a mass term in the Dirac equation, leads to 
a breaking of the TRS in the absence of the magnetic flux. 

In this paper, motivated from the above result, we study the Aharonov-Bohm- 
Coulomb (ABC) problem [T3 [2] in a graphene ring to investigate the effect of a 
Coulomb type potential in the form of £/r on the splitting of valley degeneracy, [f] 
Within the framework of field theory there are many possible ways in introducing 
the potential to the Dirac equation. The only criterion for a consistent coupling 
of the potential to the Dirac equation is the fact that the larger component of the 
Dirac field ip should satisfy the Schrodinger equation [p 2 /2M + £/r] ip = £tp in the 
non-relativistic Galilean limit [15] . In fact, there are at least two ways such that 
the same consistent Galilean limit is fulfilled. In the present work we consider the 
following two possibilities |16j : one is the scalar coupling and the other is the vector 
coupling. As we shall describe explicitly in next section, in the former case the £/r 
potential enters the Dirac equation as a mass term, while it enters the equation as 
an energy term in the latter case. From the viewpoint of time-reversal symmetry the 
scalar coupling is expected to break the TRS, but the vector coupling is not. In the 
following, we explicitly show that the scalar coupling indeed leads to the splitting of 
valley degeneracy, while the vector coupling does not. What is remarkable and new 
in our result is that, besides the splitting of the single- valley degeneracy, the scalar 
coupling produces a separation between the energy spectrum of one valley and that 
of the other valley, and the separation increases with the interaction strength £. 

The paper is organized as follows. In Sec. [2j starting with a brief recapitulation 
of the time-reversal symmetries in graphene, we give a qualitative argument how the 
TRS is broken by the £/r potential in the AB ring of graphene, whereby the splitting 
of valley degeneracy is produced. In Sec. [3] we solve the 2D Dirac equation for the 
scalar coupling of the potential in a graphene ring. We derive an analytical expression 
of the energy spectrum in terms of valley parameter (denoted by r) and the interaction 
strength £. These two parameters interplay to separate the whole energy spectrum of 

X Thus it includes both the electric Coulomb potential and a relativistic scalar potential. 
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one valley from the other. We then compute a persistent current in the ring. Owing 
to the separation of energy spectra the persistent current is asymmetric about the 
zero magnetic flux, which represents an essentially single-valley characteristic in the 
case of scalar coupling, known as the valley polarization [TTlUHj- In Sec. [4| we present 
the energy spectrum for the vector coupling. It is shown that, contrary to the scalar 
coupling case, the potential via vector coupling to the Dirac equation alone can break 
neither the intravalley nor the intervalley degeneracies because of the conservation of 
the TRS. We find that the potential in the vector coupling serves effectively as an 
additive constant in the energy spectrum and is completely decoupled from the AB 
effect. Finally, there will be a conclusion in Sec. [5] 




Figure 1. The Honeycomb lattice (left) consists of two sublattices, one with A 
atoms and the other with B atoms. The right is the corresponding Brillouin zone 
and its symmetry points. Note that the two symmetry points, K and K' , are 
related by K' = —K. The energy spectra near these points for conduction band 
are given by E = hvp\k ± K\, which are the conic dispersion relations, called 
valleys. 



2. TRS breaking and splitting of valley degeneracy 

In this section we first briefly recapitulate the discrete symmetry of the time reversal 
in graphene, then discuss the symmetry breaking in the ABC problem of graphene 
ring. As illustrated in Fig. [I] there are two inequivalent symmetry points K and K' 
in the BZ of a graphene, and to each point two sublattices A and B are associated. 
The Hamiltonian of a Dirac electron in graphene is conveniently described by two 
pseudospins related to the two valleys at K and K' and the two sublattices. In the 
following we use the Pauli matrices r = (t"i,T2,7>j) and a = (<7i, a 2 , 03) to denote the 
valley and sublattice degree of freedom, respectively and also use the 2x2 unit matrices 
To and (To- Based on the geometry in Fig. [T] the Hamiltonian and corresponding four- 
component spinor can be described as 

H a =T a ®H K , H K = V F (T-p i * = ($m,*m,-*M',*M') T (1) 

where vp ~ c/300 is the Fermi velocity, T stands for transpose, and we have used the 
valley-isotropic form of hamiltonian for convenience of subsequent calculations [17] . 
To see the TRS in graphene we introduce the following time-reversal operator 

T =[{it 2 )®{i ( j 2 )]C 1 (2) 
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where C is the complex conjugate operator. The effect of T on the Hamiltonian Hq 
and the state ^> are 

TH T~ l = H Q , T* = ^* AK ,^* BK „-** BK ^* AK ) T . (3) 

This shows that the Hamiltonian Hq is invariant under the transformation by T but it 
interchanges the valleys; there exists a degeneracy between the two valleys (henceforth, 
intervalley degeneracy) . According to Ref. [TH] , the Hamiltonian satisfies another TRS 
under the transformation by the operator 

= r <g) CT , 0^ = ia 2 C. (4) 

With this operator one can verify 

Sffoe-^flb, &* = (**BK,-**AK,*AK>>**BK>) T - ( 5 ) 

This operator exchanges the sublattices within a single valley, but does not interchange 
the valleys as we can see in the right equation. Since the Hamiltonian Hq is 
invariant this leads to a degeneracy within a single valley (henceforth, intravalley 
degeneracy) [20]. More specifically, the operator CT transforms Hr in Eq. Q to 
©aHfc ©a = that is, is invariant as the intravalley degeneracy implies: & a 
effectively changes the sign of a and p. 

Having introduced the TRS's of graphene we now discuss the TRS breaking of 
the ABC problem in a graphene ring. For this we first consider an AB type ring of 
a graphene discussed in Ref. [9]: inner and outer radii of the ring are a + w/2 and 
a — w/2, so that the ring width is w. Since the Hamiltonian Hq has a valley-isotropic 
form it is more convenient to use Hx with two-component spinor and insert the valley 
index in an appropriate place. For the AB ring the Hamiltonian then reads [§] 

H{A) = H K (A) + TV(r)a 3 , H K (A) = v F a ■ (p + eA) (6) 

where the circularly symmetric potential rV{r)a 3 is introduced to confine a Dirac 
electron on the graphene ring: V(r) = when \r — a\ < w/2 and V(r) —¥ oo when 
|r — a\ > w/2. The index t = ±1 denotes the valleys: r = +1 for the K valley and 
t = — 1 for the K' valley. We note here that this potential is proportional to 03. 

Obviously the magnetic field breaks the effective TRS of 0^: & a HK (A)©^ 1 = 
Hk(— A). In fact, the TRS of the Hamiltonian with four-component spinor is also 
broken: T[t <8> iJ^(A)]T _1 = t <£> Hk{~A). Thus, the presence of magnetic field 
will lift both of the intervalley and intravalley degeneracies. What is interesting in 
Eq. (JgJ) is the TRS breaking by the confinement potential term when A = 0. In the 
absence of magnetic field one can see 

e*[H K (0) + rVir)^}®- 1 = H K (0) - rV(r)a 3 . (7) 

This implies the intervalley degeneracy can be broken by the confinement potential as 
demonstrated in Ref. [9]. According to Berry [21] this potential enters the 2-D Dirac 
equation by the replacement of Mc 2 — > Mc 2 + V(r)a 3 with M — for a massless 
particle. In this sense the confinement potential is regarded as a mass term in the 
Dirac equation, and it can be conjectured that a mass term proportional to 173 leads 
to a breaking of the effective TRS in an AB graphene ring. 

We now turn to the problem of ABC in a graphene ring. Here we need to include 
a Coulomb type potential £/r in the Dirac equation. To address the TRS breaking 

§ Throughout this paper we use the convention h = c = 1 unless otherwise specified. 
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by this potential we begin with the general form of a 2D relativistic Dirac equation 
minimally coupled to AB potential A, which can be written as 

[c/3 7 • IT + PMc 2 } ip = Eip, (8) 

where M is the bare mass of a Dirac electron, ip is the two-component spinor, and 
Ilj = —idi — eAi is the covariant derivative multiplied by —i. The Dirac matrix is 
chosen as 

P = *3, fai = (<ri,8*2), { 7 ^,7 I/ }=2<^, if v = (+,-,-), (9) 

where s is twice of the spin value (+1 for spin "up" and —1 for spin "down"). We 
employ a thin flux tube as the AB potential in the form 
ocean 



eAi = 



(10) 



which yields the magnetic field B = —2-Ka5(r)/e along the ^-direction as expected. 
Therefore, the parameter a represents a magnetic flux in unit of —e/(2ir). 

As mentioned in introduction there are at least two possibilities for the potential 
£/r[j[|to be included in the Dirac equation, the scalar coupling and the vector coupling. 
For the scalar coupling it couples to the equation as 

r , 



ip = Eip, 



c^7 ■ n + /3 [NIc 2 
and for the vector coupling the equation becomes 
[c/?7 • II + /3Mc 2 l tp = ( E — - ) ip. 



(11) 



(12) 



From a physical point of view the potential f/r in the vector coupling corresponds 
to the electric Coulomb potential (time component of the relativistic four vector) and 
the potential in the scalar coupling is a relativistic scalar potential (other than the 
electric Coulomb potential). For the Dirac electron in a graphene, since M = and 
c = vf, we have 



H v 1> = Ei>, 



H S (A)=H K (A)- 
H V (A)=H K {A) 



-0-3, 



r 
r 



(13) 
(14) 



where H S (A) and H V (A) stand for Hamiltonians of scalar coupling and vector 
coupling, respectively, Hjc (A.) is given in Eq. and we have used the conventions 
in Eq. (^f with s = +1. 

As before the magnetic field will break both of the TRS's for T and 0. To 
see the effect of the potential £/r, we take the time-reversal transformation on the 
Hamiltonians for A = 0. For the operator T we get 



Tto^Hs^T- 1 
Tt q ®H v {<S)T- 1 



T <g) 



TO 1 



H K (0) 
H k (0) 



-0-3 



(15) 



|| It should be noted that this is a three-dimensional expression of the Coulomb potential. The true 
2D expression of the Coulomb potential is proportional to lnr. The reason for use of £/r is due to 
the fact that the two-dimensional graphene is embedded in the three-dimensional space. 
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The scalar coupling breaks the TRS for T, but the vector coupling preserves the TRS. 
For the effective TRS operator & a the two Hamiltonians are transformed to 

©^(O)©; 1 = H K (Q) - i<r 3 , 

" (16) 

&M0)®- 1 =H K (0) + ±. 

r 

Here we have the same results: the scalar coupling breaks the effective TRS for © ff , 
but the vector coupling does not. Therefore we expect that the scalar coupling will 
lift both of the intervalley and intravalley degeneracies. On the other hand, the vector 
coupling breaks neither of the two degeneracies. In the following sections we will 
compute the energy spectra of the ABC problem in a graphene ring to show the 
effects of each coupling explicitly. 



3. ABC problem with scalar coupling 



In this section we compute the energy spectrum for the scalar coupling of ABC in a 
graphene ring. As described before the geometry of the ring is the same as discussed 
in Ref. [9]. Using the Hamiltonian (13) and the confinement potential in Eq. (JgJ) the 
2D Dirac equation is 



H s 4> = Ei/>, H s 



vf<t ■ (p + eA) 



03 



rV(r)a 3 , 



(17) 



tp = r(n ± • a)ip, 



where r — ± is the valley index. As defined earlier, in the confinement potential, 
V(r) = when |r — a\ < w/2 and V(r) —> oo when |r — a| > w/2. According to 
Refs. Tfj and [22] the boundary conditions on the two-component spinor ip can be 
expressed as the form 

(— sin#, cos6>), r — a + ^ , . 

-(- sin 6, cos 9), r = a-f 1 j 

where nj_ is the unit vector perpendicular to the normal direction n on the boundary 
(i.e, n± ■ n = 0) and lies in (x, y) plane, that is, on the ring plane. This choice of 
boundary conditions is based on the requirement of the hermiticity of the Hamiltonian 
H within the boundary, which yields the condition of no outward current at any point 
on the ring boundary (i.e., n ■ v = 0, where v =< ip\vF&\ip >)-|5] 

For a Dirac electron inside the ring the Dirac equation reads, using V(r) = 0, 



cr ■ (p + eA) + -er 3 
r 



ip = Eip 



(19) 



1 More precisely, the boundary conditions are determined by the self-adjointness of the Hamiltonian 
> — < Hip\-ip >= 0. Mathematically, this is related to the self-adjoint extension 
with deficiency indcies (2, 2), so that the two-component spinor satisfies tp = Utp at boundaries, where 
U is a 2 X 2 unitary, hermitian matrix with unit determinant 23 24 . f n the present case U = n± ■ cr 
with which the operator & a = iaiC anticommutes, that is, {U, © CT } = 0. Thus, the chosen boundary 
condition does not preserve the effective TRS. This particular choice of U is to prevent the Klein 
tunneling. 
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where E = E/vf and £ = Operating [— er- (p + eA)+a^/r+E] on the equation 

and using polar coordinates we have second order equations 



1 



r 



(dg + ia f 



E 2 + seBa 3 



4> 



-isB 



(20) 







where B and a are the magnetic field and flux given in Eq. ( 10 ). Since the Hamiltonian 
(17 1 satisfies [H, J3] = 0, where J3 = — i0e + {s/2)a 3 , the solution to the Dirac equation 
can be written in the form 



1 3 

± 2' ± 2'' 



(21) 



Inserting this into Eq. (p0| we can extract the radial equation 



1 



Or 



(m + a) 2 +£ 2 + 1/4 



X2m(r) 



1 



O703 - £0-2) 



Xim(r) 
X2m(r) 



where 



r\ = m + a. 

Using the matrix diagonalization the right hand side can be written as 

<j) 4> 

U c = cos —(To + 1 sin — cri 
where e(x) — \x\/x is the alternating function, do is the 2x2 unit matrix, and 

4> . . Jv 2 + £ 2 - \v\ 

tan - = e{r]) — 




(22) 



(23) 



(24) 



(25) 



With this diagonalization the radial equation ( 22 1 reduces to 

(m + a) 2 + £ 2 + 1/4 



Of. 



1 



Or — 



E' 



f2m(r) 



1 



/2m (r) 



(26) 



where the two components fi m (r) and /2 m (r) are related to the spinor x(r) as 



X(r) 



Xlrn(r) 

X2m(r) 



u! 



fim(r) 
hm{r) 



(27) 



The reduced equation for each component in Eq. ( 26 1 is the Bessel's equation and, 
introducing dimensionless radial variable p = \E\r, the solutions can be expressed as 

hm(p) = A lm Hi l }{p) + B lm Hl 2 ](p) 

h m {p) = A 2m Hi 1 }{p) + B 2m H^(p), 
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where H^\p) and i^ 2) (p) are the Hankel functions, and the orders are given by 



L 2 + V 2 ± 1 _ e{v y 



(29) 



Substituting these solutions with the relation (27) into the spinor ip in (21 1 and using 
the Dirac equation ( 19 ) together with the recurrence relations of the Bessel equations 

(30) 



one can also derive the following relations between coefficients 
-42m = ie(r]E)A lm , B 2m = ie(rjE)B lm . 



The cigenspinor x{p) f° r the radial equation (22) is then obtained to be 



Xlm(p) 



A lm Hi 1 }{p) + B ln M 2 }{p) \ 
ie(r,E)A lm Hl 1 Hp)+ l e(r 1 E)B lm Hi 2 J(p) { > 



To determine the coefficients Ai m and B\ m for each m we use the boundary 
conditions given in Eq. (18). For the components of the spinor solution (21) the 
boundary conditions give 



X2m(Pi/\E\) = -irxim(pi/\E\), 
X2m{P2/\E\) = irxi m (P2/\E\), 
where pi = \E\(a — w/2) and p 2 = \E\(a + w/2). Defining 

Y^(v u v 2 -p) = tan - e(r,E)H^(p) 

Z W K vi;p) = H${p) + e( V E) tan |ff«(p) 
with i = 1 or 2, the boundary conditions (32) read 



(32) 



(33) 



Blrn — 



{v- ,v+;p2) + tZ (y \v_ , v+ ■ p 2 )J Ai r , 
+ [Y^{v_, V+ -p 2 )+TZW{^,u + -p 2 ) 
yW(i/_,i/+;pi)-TZW(i/_,i/ + ;p 1 )] A lm 
+ [yW(v-,v + ; Pi )-tZW(v-,u+; Pi )\ B lm = 0. 
The secular equation requires then the following relation 
Y « ,v+;P2)+ tZW (v- ,v+;p 2 ) Y^ {v_ ,u+ ; p 2 ) + rZ^(y_ ,u + ;p 2 ) 



(34) 



(35) 



v+; pi) - tZW(u-,u+; Pl ) YW(v_,v + ; Pl ) - rZW(i/_, Pl ) ' 

From this the spectrum of energy eigenvalues of a Dirac electron in the graphene ring 
can be calculated. It should be noted here tha t the scalar coupling is implied in the 
orders v± of the Hankel functions (see Eq. (29)) and hence the Eq. (35) with £ = is 
identical to the energy eigenvalue equation derived in Ref. [9]. 

To obtain an explicit expression of the eigenvalue spectrum we assume 

w V F ,„„s 

2a \E\a 

and exploit the asymptotic formula of the Hankel functions for large p in the condition 
( 35 1 . This gives the following equation 



\E\=e n + V -^e(r,)^if+e 
w ' 



VF 

\E\a 



^T, m ,s(a;^;E) 



(37) 
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where e„ = vp(n + l/2)/w (n = 0, 1, 2, . . .), and 



w 2 = 



r/ 2 + e - e(E) 



VF 

2\E\a 



(38) 




T€(E) 



Solving Eq. (37) by iteration and keeping leading terms, the energy eigenvalues are 
obtained to be 



E nm (T,^,a) = ±e n ± £ nm (r,f, a), 



where 
£nm(T,£, a) 



2a 



w 



1 

2Y r , 



(?)' 



<v)Vv 2 + £ 2 O r , m , s (a;£;±e n ). 



(39) 



(40) 



In Fig. [2] we plot the first few positive energy levels with n = as a function of 
the magnetic flux a for different values of £. In the figure the solid and dotted lines 
correspond to r = 1 and r = — 1, respectively. Fig. [2]ja) is the plot of a-dependence 
when the Coulomb interaction is zero, which is the same situation as Ref. [§]■ For the 
analysis of the energy spectrum with £ = 0, we let E O m(T,0,a) = E(r,m,a). When 
a = 0, since the time reversal symmetry by T is preserved, there are degeneracies 
between E(l, m, 0) and E{— 1, — m, 0). However, the intravalley degeneracy is broken, 
that is, E(t, m, 0) ^ E(t, — m, 0) because of the confinement potential that breaks the 
effective time reversal symmetry by 0. Obviously, when q ^ 0, these degeneracies 
are broken due to the AB potential. 

Fig.^b) and Fig.[2]jc) are the plots of a-dependence when £ = 0.2 and £ = 0.4, the 
case of repulsive interaction. As noted from the figures, the splitting of the intravalley 
degeneracy still exists due to the mass terms of the scalar coupling and the confinement 
potential. What is more interesting here is that the scalar coupling raises the whole 
energy spectrum of the r = +1 valley and lowers the whole energy spectrum of the 
r = — 1 valley. The situation will be reversed if the interaction is attractive; the 
t = — 1 valley is raised, while the r = +1 valley is lowered. As explained in Sec. [2] 
this effect is ascribed to the TRS breaking by the scalar coupling of £/r entered as a 
mass term in the Dirac equation. 

Another interesting feature to be noted is that the separation between the two 
spectra gets larger as the interaction strength increases. To see the energy difference 
between the two separated spectra, we look into the r-dependent part in the Eq. ( 40 1 



for a fixed set of (n,m,a). The energy difference is then determined by the term 
^(t^q) which can be positive or negative depending on the sign of r together 
with r\ and uj\. With the assumption (36), since rjcdi ~ ($ a positive constant), 
we may write £ nm (r,£,a) oc e(T£)£y r\ 2 + £ 2 , which becomes positive when r£ > 
and negative when t£ < 0. Thus, for a repulsive interaction for which £ > 0, the 
energy spectrum of r = +1 valley is raised, while that of r = — 1 valley is lowered; 
for the attractive interaction, since £ < 0, the opposite occurs. The energy difference 
between the two spectra is then A£ nm oc rf + £ 2 , which explains the increase of 
the separation with £. 

A possible way of observing the splittings of valley degeneracies is to measure 
the magnetic moment induced by persistent current in the graphene ring |26) . 
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5 = 






m=5/2 / ~ 


^\m=-3/2 _ . ' 






m=3/2^> < C. 











(b) 



1.58 
1.57 
1.56 



. ? = 0.2 










Figure 2. Energy spectrum of positive Eo m (r, £, a). For convenience we choose 
Vf = w = 1 and a = 10. The solid and dotted lines correspond to t = 1 
and t = — 1, respectively. The colors pink, orange, red, green, blue, and cyan 
correspond to m = 5/2, 3/2, 1/2, —1/2, —3/2, and —5/2, respectively, (a) is the 
energy spectrum when £ = 0, the case without interaction which is identical to 
Ref. [9]. (b) and (c) are the energy spectra when the interaction exists. Note that, 
besides the breaking of valley degeneracy, the spectrum of the r = +1 valley is 
raised and the spectrum of the r = — 1 valley is lowered. 



The persistent current can be obtained from the energy eigenvalues and, at zero 
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- 5 = 




N = 4 








N = 3 






1 N = 2 




N = 1 





(b) 



(C) 12 




Figure 3. Persistent currents as a function of a for different values of £. N is the 
number of electrons in the occupied states, (a) is when 5 = 0, the case without 
Coulomb interaction, (b) and (c) are the cases with Coulomb interaction through 
the scalar coupling. Note that persistent currents are symmetric about a = 
when 5 = 0, but asymmetric about a = when 5 > 0. The latter is ascribed to 
the separation between the spectrum of r = +1 valley and that of r = — 1 valley. 



temperature, it is given by 

J(€,a) = - EE^>"»( T >^)- (41) 

t=±1 n rn 
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In Fig. [3] we present the persistent current as a function of the magnetic flux a for 
different values of £ and the number of electrons N (including spin). Fig. [3](a) is the 
plot when £ = 0, the case without interaction. Since the intervalley degeneracy exists 
the current is symmetric about a — when there are equal number of electrons at 
each valley, as can be seen from the figure. When £ > and large, the separation 
between the two spectra of r = ±1 valleys is large. Since the electrons occupy from 
the lowest levels at zero temperature the states of the lower spectrum (the r = — 1 
valley) will be occupied first, while those of the higher spectrum (the r = +1 valley) 
are almost empty. Thus, the persistent current is contributed largely by the lower 
valley electrons, and becomes essentially a single valley phenomenon, known as valley 
polarization. This should produce an asymmetric persistent current about a = and 
a finite value at a — 0, which are illustrated in Figs.[3jb) and[3](c). As we can notice in 
Fig.[3]Jb), even for very small value of the interaction strength, the qualitative feature 
is quite different from the case without the interaction because of its role of symmetry 
breaking. We also emphasize that the a-dependences of persistent currents for £ > 
are essentially the same because the electrons occupy only the energy levels in the 
lower spectrum corresponding to the r = — 1 valley. 



4. ABC problem with vector coupling 



In this section we consider the vector coupling of the potential £/r. Here we will 
use the term Coulomb potential because the vector coupling is related to the electric 
Coulomb potential. Using the Hamitonian given in Eq. (14) the Dirac equation inside 
the ring can be written as 



er ■ (p + eA) ip= IE 



(42) 



Operating <x ■ (jp + eA) on the equation and following the same procedure as in the 
previous section the reduced radial equation is given by 



1 



i] 2 + 1/4 



~\ 2 



E 2 



fim{r) 

/2m(f) 



(43) 



= 2 e WV'' 2 "^ (T 3 



hm(r) 
hm{r) 



Here, the two components {fi m ( r ),f2m( r )) are related to the spinor %(r) as follows: 



X(r) 



Xim(r) 
X2m(r) 



U" 1 



hm{r) 

/2mW 



(44) 



where the transformation matrix U v is defined by 



U^ 1 = cosh £ cr ± e(O e ( ? 7) smh tt °2 



tanh- = e(f) 



V 2 - £ 2 



(45) 



and satisfies the matrix diagonalization 77173 — i^a\ = ?7 2 ~ (, 2 U V 1 o~3U v . Note 

that U v in Eq. (44) is not unitary. 
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The solutions of Eq. ( 43 ) are expressed in terms of the confluent hypergeomctric 
functions: 

fi m (p) = e v (-2zp) t - [C lm M(a-, 6-, -2ip) + D lm U(a_, 6_, -2ip)} 

(46) 

/2m (p) = e lp {-2ip) t + [C 2m M(a+, b+, -2ip) + D 2m U{a+, b+, -2ip)} 

where M(a, &, p) and U{a, 6, p) are the Rummer's functions [27) . p = \E\r, and a±, b± 
and i± are given by 

a± = t± + - + ie(E)l b± = 2t± + 1 



t ± =J^f± e M. (47) 

By the same method as in Sec. II and exploiting the recurrence relations and 
differential properties of the confluent hypergeometric functions, one can derive the 
following relations between coefficients after some algebra; 

C 2m = e{E)c 21 C lm D 2m = e(E)d 21 D lm , (48) 

where 

e(»7) / ; — _ N t(v) 



V 2 -e + HE)i 



c 2i = - I „ - , r- I dai= I J ^ I ■ (49) 



^2 _£2} +2 J r) 2 -£2 



Substituting the solutions ( 46 1 with the relation ( 48 ) into the transformation equation 
(44 1 we can write the eigenspinor x(p) as 

Xim(p) ) =e i P{ _ 2ip) t- ( mr(p) uM \ / C lm \ 
X2m[p) ) \ m 2 {p) u 2 (p) ) \ D lm ) 

where 

mi(p) = cosh ^ M(o_, -2ip) + ie(£r)E)(-2ip) e( -^ c 21 sinh ^M(a + , b + , -2ip) 

ui(p) = cosh - U(a-, b-, -2ip) + ie{iriE){-2ipY {r)) d 21 sinh -U(a + ,b + ,-2ip) 

6 6 (51) 

m 2 {p) = e(E)(-2ip) e ™ C21 cos h | M(a +! 6 + , -2tp) - ie(|»/) sinh |M(o_, &_, -2*/») 

u 2 (p) = e(S)(-2«p) e (") dji cosh | C/(a+, 6+, -2ip) - ie(|r?) sinh ^(o_, -2ip). 

To proceed we use the same infinite mass boundary condition introduced in Sec. II. 
Inserting Eq. ( 50 1 into the boundary condition ( 32 ) it is straightforward to show 

[m 2 (p 2 ) - M-mi(p 2 )] Cim + [vaiPi) ~ iTUi(p 2 )] D lm = 
[m 2 (pi) + irmiipi)] C lm + [u 2 (pi) + mti(pi)] D lm = 0, 

where pi = |£7| (a — u>/2) and p 2 = (a + w/2). The secular equation requires then 
the following condition: 

m 2 (p 2 ) ~ iTini(p 2 ) u 2 (p 2 ) - irui(p 2 ) 



m 2 (pi) +«rmi(pi) it 2 (pi) + m*i(pi) ' 



(53) 
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(b) 




Figure 4. Energy spectrum for the vector coupling. For convenience we choose 
vp = w = I and a = 10. The solid and dotted lines correspond to t = 1 
and t = —1, respectively. The colors pink, orange, red, green, blue, and cyan 
correspond to m = 5/2, 3/2, 1/2, —1/2, —3/2, and —5/2, respectively, (a) and 
(b) are energy spectra for £ = 0.2 and £ = 0.6, respectively. Note that the 
intervalley symmetry about a = is still kept here even in the presence of the 
potential and the whole spectrum is raised by £/a compared to the Fig.[2ja). 



To obtain the eigenvalue spectrum we impose the same condition with Eq. ( 36 1 and 
use the asymptotic formula £J of the confluent hypergeometric functions for large p. 
Expanding Eq. ( 53 ) up to the second order of perturbation, one can show that the 
energy eigenvalue satisfies 



\E\ = 



10 



VF 

\E\a 



, (54) 



where n is a non-negative integer. Solving Eq. (54 1 by iteration and keeping only 
leading terms, we obtain the energy eigenvalues 



E nm (T, £, a) = ±e n ± X n r] r\ =F 



i 

i 

a 



(55) 



+ For large \z\ 
M(a, b, z) ~ c :" 



r(o) 



E 



(b - a) n (l - a)r, 



rr, , \ -a ( a )nl + a - b )n , , _„ 



where (a) n = a(a + 1) ■ ■ ■ (a + n — 1); see Ref. [27 
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where e n — vp(n + 1/2) /w (n = 0,1,2,...) and A„ = (vp / 'a) 2 /2e n , 
dependence is through the variable ?y defined in Eq. ( 23 1 . 



and the a 



We can immediately see from Eq. (55 1 that the effect of the potential £/r in the 



vector coupling is merely an additive constant (£/a), shifting the whole spectrum of 
energy eigenvalues. It also reveals the complete decoupling between the AB and the 
Coulomb effects, because the Coulomb effect appears in the first order whereas the AB 
effect arises in the second order perturbation. In Fig. |4j we show the energy spectrum 
in which the symmetry between the two valley spectra about a — is manifested, 
representing the existence of the intervalley degeneracy. As discussed in Sec. [2] the 
Hamiltonian without AB potential in Eq. ( 42 1 is invariant under the time reversal 



operators T and 0, and hence the vector coupling breaks neither the intervalley nor 
the intravalley degeneracies. It is thus expected that the eigenvalue spectrum in the 
vector coupling has essentially the same feature as the spectrum without the Coulomb 
interaction. By the same reason, the Coulomb potential in the vector coupling does 
not alter the behavior of the persistent current and hence it is symmetric about a = 
as shown in Fig. ^&) . 



^ scalar coupling 

















-0.4 -0.2 0.0 0.2 0.4 



(b) 



vector coupling 





























-0.4 -0.2 0.0 0.2 0.4 



Figure 5. The ^-dependence of the energy spectrum when there is no AB 
potential (that is, when a = 0). For convenience we let Eo m (T, £, 0) = E(m,T,£), 
and the parenthesis denotes (m, t). (a) and (b) correspond to the scalar coupling 
and the vector coupling for m = ±1/2, ±5/2 and r = ±1, respectively. The 
splitting of intervalley degeneracy is clearly seen in (a), but it is not broken in 
(b). 
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In Fig. [5j we also present the energy spectrum as a function of £ without the 
AB potential to compare the bare effects of the Coulomb type potential between the 
scalar coupling and the vector coupling. Fig. [5^ a) shows the splitting of intervalley 
degeneracy due to the symmetry-breaking term in the scalar coupling. However, in 
Fig. [5^b), we see that no intervalley splitting exists in the case of vector coupling. 
This also confirms that, in the case of vector coupling, the Coulomb potential alone 
does not break the intervalley degeneracy, as explained by the consideration of the 
time reversal symmetry. 

5. Conclusion 

In this paper we have considered the ABC problem in a graphene ring with magnetic 
flux tube and a Coulomb type potential £/r at its center. We have investigated 
the effects of the potential on the energy spectrum in two different ways: the scalar 
coupling and the vector coupling. In the scalar coupling the potential enters the 2D 
Dirac equation as a symmetry-breaking mass term, so that both the intervalley and 
the intravalley degeneracies are broken. The main effect of the interaction appears as 
the separation between the energy spectrum of one valley and that of the other valley, 
which is attributed to the breaking of the intervalley degeneracy; for the repulsive 
interaction the t = +1 (K) valley is lifted while the r = — 1 (K') valley is lowered, 
and vice versa for the attractive interaction. Contrary to the scalar coupling, the 
potential in the vector coupling does not break any symmetry and only shifts the 
AB energy spectrum, indicating essentially the same feature as the energy spectrum 
without a Coulomb potential. 

The results obtained here suggests that the scalar coupling of the Coulomb 
type potential £/r can decouple the two valleys because of the considerable lift of 
energy spectrum, so that each valley degree of freedom becomes an independent 
quantity, called a valley isospin. A real experiment with graphene ring may realize 
this decoupling, known as the valley polarization, which is an important element in 
the graphene-based electronics, called valleytronics [T51 128] . 

Another interesting issue related to the AB effect is when the thin magnetic flux 
tube is located at the center of the honeycomb lattice shown in Fig. [lja). Since the 
anyon is a flux-carrying particle, a similar situation arises when there is an anyon 
impurity in the graphene plane. Since, in this case, the magnetic flux is located on 
the graphene, we should treat the singular solution problem, which was extensively 
discussed in Refs. [5^1 [SO] in the context of the anyonic and cosmic string theories. 
Roughly speaking, there are two prescriptions for the interpretation of the singular 
solution: one is a mathematically-based prescription called self-adjoint extension and 
the other is a physically-based prescription. We do not know yet which prescription 
is physically more reasonable. Probably, the real experiment with graphene may shed 
light on the treatment of the singular solution. If so, our understanding on the anyonic 
and cosmic string theories can be greatly enhanced through the graphene physics. We 
will explore this issue in the future. 
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